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ABSTRACT 

Asteroseismology of F-type stars has been hindered by an ambiguity in identification of their oscillation 
modes. The regular mode pattern that makes this task trivial in cooler stars is masked by increased 
linewidths. The absolute mode frequencies, encapsulated in the asteroseismic variable e, can help solve 
this impasse because the values of e implied by the two possible mode identifications are distinct. We 
find that the correct e can be deduced from the effective temperature and the linewidths and we apply 
these methods to a sample of solar-like oscillators observed with Kepler. 
Subject headings: stars: fundamental parameters — stars: interiors — stars: oscillations 



1. INTRODUCTION 

Asteroseismology of solar-like stars is developing 
rapidly, driven by the succes ses of the spac e telescopes 
CoRoT (IMichel et all [2001 and Kevler dKoch et all 
[20101 IGilliland et all 120101: IChaplin et all 1201 ID . By 
studying the oscillation modes of these stars, infer- 
ences can be mad e about their interior structures (e.g. 
IVerner et al.ll20iT[) . Except for the most basic analyses, 
it is crucial to identify the oscillation modes, that is, the 
radial order n and the spherical degree I. 

In the Sun and similar stars, mode identification is 
straightforward thanks to the distinctive pattern of alter- 
nating odd and even modes in the power spectrum. This 
pattern consists of a regular sequence of I = 1 modes, 
interspersed with close pairs of I = and 2 modes that 
fall approximately halfway in between. However, stars 
significantly hotter than the Sun have large linewidths 
that blur the I = 0, 2 pairs and make mode identifi- 
cation very difficult. In this Letter we demonstrate a 
solution t o this problem that applies the method pro- 
posed by iBedding fc Kieldsenl (|2010D and iWhite et al] 
(|2011af) , which uses the absolute frequencies of the os- 
cillation modes. 

2. METHODS 

For main-sequence stars, the frequenci es are well- 
approximated by the asymptoti c relation (jVandakurovl 
IT967t lTass O1 illT980l : iGoTTghlllQSfll 

v nt i « Al> (n+ ~ +e) - 6v i- (1) 



Here, Av is the large separation between modes of the 
same degree I and consecutive order n, Svqi is the small 
separation between modes of different degree and e is 
a dimensionless offset, which we discuss in greater de- 
tail below. Typically, only modes of I < 2 are ob- 
served in intensity due to cancellation over the unre- 
solved stellar disk, although I = 3 modes can be ob- 
served in the highest signal-to-noise targets, such as 11 
Kepler subgi ants for which frequencies have been de- 
termined by lAppourchaux et all (I2012bf) and the solar 
analogs 16 Cvg A and B iMetcalfe et all (|2012ft . 

The asymptotic relation makes it easy to determine 
the mode degrees for the Sun and similar stars. Each 
I = mode is separated by 5vq2 from an I = 2, and 
separated by Av/2— Svqi from the I = 1 mode of the same 
order. An example is shown in Figure for the Kepler 
star KIC 6933899, which has an effective temperature of 
5840 K. 

Mode linewidth inc re ases wi th effective te mper- 
ature (IChaplin et all 120091: iBaudin et all l20TTt 
lAppourchaux et al.l l2012at iCorsaro et al.l I2012I ). re- 
flecting shorter mode lifetimes in hotter stars. In 
some F-type stars, the linewidths become so large that 
the pairs of I — and I = 2 modes are unresolved 
and it becomes difficult to distinguish between the 
blended I = 0, 2 modes and the I = 1 modes. This 
problem w as first observed by CoRoT , in the F5 star 
HD 49933 (|Appourchaux et all I2008D and has since 
been seen in othe r CoRoT stars (jBarban et al.l 120091 : 
iGarcia et al.l I2009D an d in the bright F5 star Procyon 
(|Bedding et al.l I2010H) . We also see it in many Kepler 
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Fig. 1. — Power spectra of (a) a G star, KIC 6933899, and (b) an F star, KIC 2837475, with their corresponding echcllc diagrams (c) and 
(d), respectively. The red curves show the power spectra after smoothing. Mode identification of the G star is trivial, with modes of I = 
(orange), 1 (blue) and 2 (green) labeled. For the F star it is not clear whether the peaks labeled 'A' (blue) or 'B' (orange) correspond to 
the I = 1 or I = 0, 2 modes. 



stars and Figure [T}d shows one example, KIC 2837475 
(T eff = 6690K). 

One way to resolve this identification problem 
is to fit both possible mode identifications and 
compare the relative likeli ho ods of the two sce- 
narios (lAppourchaux et al.l 120081: IBenomar et al .l 
2001 iGruberbauer et al.l l2009t iKallinger et all l201( 



Bedding et al.l l2010bl: lHandberg fc Campantd l20TTfT . 
This relies on the profile of the even-l modes being 
significantly broader and also asymmetric, relative to 
the I = 1 modes (owing to the presence of the smaller 
amplitude I = 2 modes at a slightly lower frequency than 
the I = modes). The correct scenario should provide a 
better fit to the power spectrum. Difficulties arise at low 
signal-to-noise and with short observations, for which 
the Lorentzian mode profiles are not well resolved. Thi s 
method was first applied bv lAppourchaux et al.l ([2008), 
who fitted both scenarios for HD 49933. However, with 
additional dat a their preferred mod e identification was 
overturned by IBenomar et al.l (|2009T > . 

Other methods have been suggeste d that utilize the 
sign of the small separatio n Si/qi (|Roxburghl 120091 : 
iMosser fe Appourchauxl [2009T ) . In main sequence stars 
like the Sun, Svqi is known to be positive. How- 
ever, in many red giants 8vm is found to be negative 
(Bedding et all l2010al iHuber et all 120101: IMosser et all 
120111). so at som e point in the evolution the sign must 
flip ( Stellol [201 ID . To further complicate matters, the 
value of Svqi is quite small. At low signal-to-noise, it 
may be difficult to obtain frequencies precise enough to 



d etermine the sign of 6 Vqi re liably. 

iBedding fc Kieldsenl ()2010[) have suggested that scal- 
ing the frequencies of a star with a known mode identi- 
fication could reveal the correct mode identification in a 
second star. This method seeks to use information con- 
tained within the value of e. For this to be effective, 
e must vary slowly as a function of stellar parameters. 
This is indeed the case, with a tight relationship between 
e and effective temp erature, T e g, found both in models 
(White et al. 2011b) and in observations of Sun-like stars 
(Wh ite et al.ll2011af ). Thus, e promises to be an effective 
way to determine mode identifications, since the differ- 
ence in the value of e for the two possible scenarios is 
large (0.5). 

The existence of a relation between e and T e g is 
not surprising. The value of e is determined by 
the upper and lower tu rning points of the acoustic 
waves (e.g. iGoughl I1986J ). As such, e is heavily de- 
pendent upon the stellar atmosphere, of which T e g is 
a significant parameter. Due to inadequate model- 
ing of the near-surface layers, there is a well-known 
offset between obser ved and computed oscillation fre- 
quencies in the Sun flChristcnsen-Dalsga ard et al.l 1198 
Dziembowski et all 119881: iC hristcnscn-Dals gaard et al 



1996t iChristensen-Dalsgaard fc Thompson! 1997) and 
also in other stars (iKieldsen et al.1 l200"8FlWhite et al.l 
PiOllat IMathur et al.l I2012TI . This offset results in the 
computed e being smaller than observed, typically by 
~ 0.2 as inferred from the displacement of model tracks 
from observations in the e diagram. 
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Fig. 2. — Modified H-R diagram: average large frequency sep- 
aration, Av, against effective temperature for stars in our sam- 
ple. Stars with secure mode identifications are indicated by black 
circles. Those without are r ed squares. Grey lines are ASTEC 
(Christcnscn-Dalsgaard 2008) evolutionary tracks for a metallic- 
ity of Z = 0.011 ([Fe/Hl = - 0.2 de x), matching the T cS cal- 
ibration of Pinsonncault ct al. (2011). The dashed line indi- 
cates approximately the cool edge of the classical instability strip 
USaio fc Gautschy|[l99l 'l. 

The purpose of this Letter is to extend the relationship 
between e and T e g to higher temperatures, and thereby 
make reliable mode identifications in F-type stars. 

3. OBSERVATIONS AND DATA ANALYSIS 

We used observations of solar-like oscillations in 163 
stars taken with the Kepler space telescope between May 
2009 and March 2011 (Quarters 1-8). Each star was ob- 
served in Kepler's short-cadence mode (58.9 s sampling) 
for part of this period. The time series were prepared 
from t he raw observations as described by I Jenkins et al.l 
(2010) and further c orrected to remove outliers and 
jumps as described bv iGarcia et al.l (|2011ft . 

Effective temperatures were dete rmined from SPSS 
qriz c olor-temperature relations by iPinsonneault et all 
(|2011ft . Spectroscopic temperatures have also been deter - 
mined for 77 stars in our sample bv lBruntt et al.l (2012). 
In almost all cases the photometric and spectroscopic 
temperatures agree, except for several stars where the 
disagreement may be due to differing metallicities or 
unresolved binaries (the temperatures for KIC 3424541, 
3456181, 4638884, 6679371, 7976303, 8938364, 9908400, 
10018963, 10124866 and 10963065 were found to disagree 
by more than 3cr). Figure [2] shows a modified HR dia- 
gram of this sample, where we have used large separation 
instead of luminosity. 

To measure the value of e for each star, we first de- 
termined the frequencies of the I = modes. Where it 
was possible to resolve the I = and I = 2 modes, we 
measured the frequencies of the I — modes from the 
peak in the power spectrum after smoothing. This was 
possible for 115 stars. An example is shown in FigureQJt. 

Where it was not possible to resolve the / = and 
I = 2 modes, and therefore not possible to easily deter- 
mine the correct mode identification, we determined e 
for both scenarios (43 stars). In this case we used the 



frequencies of the ridge centroids, determined from the 
peaks of the heavily smoothed power spectrum, as shown 
in an example F star in Figure [T]d. There were also 
two cooler stars for which the mode identification was 
not clear (KIC 11401708 and 12555505); these have low 
signal-to- noise and the I = 2 modes are not apparent. 

For five stars with blended I — and / = 2 modes it 
was still possible to make an unambiguous mode iden- 
tification (KIC 6064910, 6766513, 7668623, 7800289 and 
8026226). In these stars, avoided crossings 'bump' the 
I = 1 modes from their asymptotically expected posi- 
tion (jOsakil [19751: lAizenman et al.lll977j ). revealing the 
correct identification. In these five cases, the frequencies 
were obtained from the centroids of the I = 0, 2 ridge. 

The values of Av and e were obtained from a weighted 
lea st-squares fit to the I = frequencies, as described 
bv lWhite et all (1201 lbl) . The weights were given by a 
Gaussian function centered at the frequency of maximum 
power, ^max, with a FWHM of 0.25 ^ ma x- By the asymp- 
totic relation, equation [T] the gradient of this fit is Av 
and the y-intercept is tAv. 

To confirm the validity of our method of frequency 
determination, we compared the values of e derived from 
our frequencies to those derived from frequencies that 
have been determined by more traditional 'peak-bagging' 
methods in the 61 stars fo r which this has been done 
(|Appourchaux et aT1l2012bft . We found good agreement, 
with e values typically agreeing to within 0.1, which is 
approximately the size of the typical uncertainty. This 
agreement is of particular importance for hotter stars 
because using ridge centroids instead of I — frequencies 
potentially biases e towards slightly lower values. We 
found this bias to be negligible within the uncertainties. 

4. RESULTS 

In Figure [3^ we show the so-called e diagram for the 
120 Kepler stars whose mode identifications were unam- 
biguous. The observed values of e are clearly offset to the 
right of the models which, as mentioned above, arises 
from the improper modeling of the near-surface layers. 
From this figure it appears that the offset in e may be 
roughly the same for all stars, which corresponds to a 
fix ed fraction of the larg e separation, as was also inferred 
bv lMathur et all (|2012l ). 

In Figure [4^ we show e versus T e ff for our sample. The 
Sun is marked in black by its usual symbol, and the stars 
with unambiguous mode identifications are colored blue. 
The trend of decreasing e with increasing T c g is clearly 
seen in these stars. 

For the 43 Kepler stars whose mode identifications are 
uncertain, the relationship between e and T c g can help. 
In Figure HJi these stars are shown in gray for both sce- 
narios. Owing to the potential ambiguity in n, each sce- 
nario is also plotted shifted by ±1. We also include bot h 
scenarios of the F stars Procyon (iBedding et al.ll2010bft . 
HP 49933 (lAppourchaux et all 12008ft and HP 181420 
(|Barban et al.H2009ft . 

For most of the stars with ambiguous mode identifica- 
tions, only one of the scenarios lies along the e-T e g trend 
defined by the stars with secure identifications and we 
adopt this as the correct one. For the pre yiously studied 
F star s, we prefer Scenario B of P rocyon (Bedding et al. 
I2010bft . Scenario B of HP 49 933 (iBenomar et aljl2009ft 
and Scenario 1 of HP 181420 (|Barban et al.ll2009f) . Pue 
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Fig. 3. — The e diagram: large separation, Au, against e. (a) Only stars with secure mode identifications are shown (filled black 
circles). Lines are ASTEC evolutionary tracks, as shown in Figure [2] although for clarity, segments of the tracks which are hotter than the 
approximate cool edge of the instability strip are not shown. Note the offset between models and observations, (b) Symbol colors reflect 
the measured effective temperature of the star. Stars with obvious identifications are circles, and those for which we can reliably make the 
identification from the e— T e g relation are indicated by diamonds. The possible identifications of one star for which the identification is still 
ambiguous in the e— T e g plane is indicated by the encircled squares. Comparing the possible identifications with the temperatures of stars 
of similar Av, the scenario on the left can be preferred. 

to the width of the e-T e ff relation, there are a few stars 
for which the situation is still somewhat ambiguous (red 
circles in Figure |4ji) . The two scenarios in these stars 
have values of e that fall towards the top and bottom 
of the relation. For these stars we must use additional 
information to resolve the ambiguity. 

To overcome the spread in the e— T e g relation, the value 
of Av is useful. For a given effective temperature, stars 
with higher masses have smaller values of Av, as can be 
seen from the models in Figure [2] Models also indicate 
that higher-mass st ars have a smaller e (see Figure^ and 
iWhite et aLll2011bf ). It follows that in the e-T eff plane, 
stars that fall towards the bottom of the trend will be 
more massive and should therefore have smaller values 
of Av than lower-mass stars of similar temperature. 

We illustrate this in the Av-e plane in Figure[3}), where 
we show stars in which the identification was already 
obvious, as well as those for which the identification could 
be readily made from the e-T e g relation. Symbols are 
colored according to T c g, with the trend of decreasing 
T e ff with increasing e quite clear. The gradual decrease in 
e with decreasing Av along lines of constant temperature 
is also visible. By comparing a star whose identification 
is still ambiguous with stars of a similar large separation, 
we can better make a decision on the identification. We 
show the example of KIC 11290197, whose two scenarios 
are circled in Figures [3}) and 0^- By comparing its T e ff 
with stars of a similar Av, we establish the preferred 
scenario, which in this case has the lower value of e. 

Other seismic parameters may also be useful. The 



method of M osser et al.l (|2010D . which uses the value of 
the small separation 8vqi, agrees for all stars, except one 
(KIC 5431016), which has a low signal-to-noise. 

Another very useful parameter is mode linewidth (T). 
As mentioned above, linewidth is strongly correlated 
with effective temperature, and so there should also be 
a correlation between linewidth and e. Linewidths for 41 
of our stars (including 12 with a mbiguous identifica- 
tions) w ere previously measured by Appourc haux et "all 
(|2012a|) . Using the SYD method described in that paper, 
we have measured linewidths in a further 26 stars with 
ambiguous identifications (excluding the five ambiguous 
stars with the lowest signal-to- noise) . Figure [4}d shows 
the relation between T and e for these 67 stars. 

We can quantify the likelihood of one scenario over 
the other by comparing how far each lies from the var- 
ious relations. To do so, we first performed a Bayesian 
linear fit to the e-T e s and e-ln(r) relations for all stars 
with T g > 5800 K with unambiguous identifications (i.e. 
the blue points in Figures 2^, and b). For all stars with 
an uncertain identification, we then calculated the like- 
lihood of obtaining the observed T e g for each of the two 
possible values of e (ea and ee), the observed T, the pa- 
rameters of the linear fits, and all respective uncertain- 
tiefl Full details of this method will be provided in a 
future paper (M. Gruberbauer et al. 2012, in prep.). The 



Gaussian uncertainties were assumed for T e g, e and ln(T); 
uncertainties for the linear parameters were determined from the 
marginal posterior of the linear fit. 



MODE IDENTIFICATION IN KEPLER F STARS 



5 




4500 5000 5500 6000 6500 7000 
T.„ (K) 




1.6 



1.4 



1.3 



1.0 



0.8 - 



0.6 - 




4500 



5000 



5500 6000 
T.« (K) 



6500 



7000 



Fig. 4. — (a) e as a function of effective temperature. Stars with 
unambiguous mode identifications are indicated by blue circles. 
Stars with ambiguous identifications have two possible values of 
e (gray circles) corresponding to the two possible identifications. 
Several stars for which the width of the relation between e and T c g 
makes the identification difficult are red circles. The scenarios of 
an example star for which we use Af to aid in the identification is 
circled, as it is in Figure |3p. (b) Same as (a), but showing e as a 
function of mode linewidth, T. (c) The selected identification for 
all stars shown in the e-T e g plane. The outlier is KIC 17258f 5, for 

which r and Tj^ 10 * disagree (see text). In all panels the Sun is indi- 
cated by its usual symbol. Identifications of Procyon are indicated 
by magenta diamonds; those of the CoRoT F stars HD 49933 and 
HD I8I420 are orange triangles. 



required integration over the parameter sp ace was car- 
ried out using MultiNest (|Feroz et al.ll2009f) . These two 
likelihood values, P(T e ff |eA,r) and P(T e ff |eB,r), were 
then used to calculate the Bayes factor (ratio of the like- 
lihoods), and hence the odds ratio and probability of 
each scenario, P(eA|7off,r) an d P((-b\T c s , T), assuming 
equal prior probability for both identifications. We per- 
formed this calculation using both photometric and spec- 
troscopic effective temperatures, where available. Our 
preferred identification is the one with the greatest prob- 
ability. We denote the value of e for the preferred iden- 
tification as e prc f, and use e a it for the alternate value. 

In Table Q] we list for each star its measured A^, the 
value of e prc f and e a it, T, both effective temperatures, 

T^ ot and T o s p° c , and the probabilities, P(e prcf |T o p ff hot , T) 
and P(e pro f |T c S g CC , r), of our preferred scenario. For 
most of the stars we find strong support for our pre- 
ferred scenario. One star (KIC 1725815), without a mea- 
sured T^ cc , is an outlier in the e~T c g plane. While the 

linewidth supports one scenario, T^ ot supports the al- 
ternate scenario. It is not entirely clear which is cor- 
rect, although the method of iMosser et al.l (|2010f ) agrees 

with the alternate scenario favored by T^g 0t . Spectro- 
scopic measurement of the temperature may help resolve 
this case. 

The e-T c g relation for the final selected identification 
of all stars in our sample is shown in Figure 2]:. 



5. CONCLUSIONS 

We have presented a method to effectively determine 
the correct mode identification in stars for which this has 
previously been a problem. These are the F stars with 
large linewidths that make it difficult to distinguish the 
I — and I — 2 modes. This method uses the relationship 
between effective temperature, mode linewidth and e to 
determine what values of e are reasonable for the star, 
and therefore which of the two possible scenarios is most 
likely correct. This method provides robust results in the 
vast majority of cases because the value of e implied by 
each scenario is very distinct, even in low signal-to-noise 
stars, representing a major improvement over previous 
methods. For the few cases that are still ambiguous, 
additional information, such as the large separation, Av, 
can be included to help resolve the matter. 
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T eS AND P(e) IN STARS WITH AMBIGUOUS MODE IDENTIFICATIONS. 
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a Identification favored by T^ ot disagrees with identification favored by T. 
b Identification made with the aid of Aia 
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